Eastern Box Turtle Movement in North Carolina

Data Collected by N.C. Wildlife Commission

Data accessed through www.movebank.org

Kernel-Density Estimation

kde_raster <- function(filename){
  data <- read.csv(file = filename)
  x <- as.data.frame(data$easting)
  y <- as.data.frame(data$northing)
  xy <- c(x,y)
  data.proj <- SpatialPointsDataFrame(xy,data, proj4string = CRS("+proj=utm +zone=17 +ellps=WGS84 +units=m +no_defs"))
  xy <- SpatialPoints(data.proj@coords)
  kde<-kernelUD(xy, h="href", kern="bivnorm", grid=100)
  ver <- getverticeshr(kde, 95)
  ver1 <- getverticeshr(kde, 75)
  ver2 <- getverticeshr(kde, 50)
  kde.points <- cbind((data.frame(data.proj@coords)),data$individual)
  colnames(kde.points) <- c("x","y","identifier")
  kde.poly <- fortify(ver, region = "id")
  kde.poly1<- fortify(ver1, region = "id")
  kde.poly2<- fortify(ver2, region ="id")
  units <- grid.text(paste(round(ver$area,2)," ha"), x=0.85,  y=0.95,
                     gp=gpar(fontface=4, col="white", cex=0.9), draw = FALSE)
  kde.plot <- autoplot(raster_utm, expand = TRUE) + theme_bw() + theme(legend.position="none") +
    theme(panel.border = element_rect(colour = "black", fill=NA, size=1)) +
    geom_polygon(data=kde.poly, aes(x=kde.poly$long, y=kde.poly$lat), alpha=0.5, color = "red", fill = "red") +
    geom_polygon(data=kde.poly1, aes(x=kde.poly1$long, y=kde.poly1$lat), alpha=0.5, color = "green", fill = "green") +
    geom_polygon(data=kde.poly2, aes(x=kde.poly2$long, y=kde.poly2$lat), alpha=0.5, color = "blue", fill = "blue") +
    geom_point(data=kde.points, aes(x=x, y=y)) +
    labs(x="Easting (m)", y="Northing (m)", title=kde.points$identifier) +
    theme(legend.position="none", plot.title = element_text(face = "bold", hjust = 0.5)) + 
    annotation_custom(units)
  kde.plot
}

pblapply(files, kde_raster)

  |                                                                                        | 0 % ~calculating  
  |==============================                                                          | 33% ~01s          
  |===========================================================                             | 67% ~01s          
  |========================================================================================| 100% elapsed=02s  
[[1]]

[[2]]

[[3]]

LS0tDQp0aXRsZTogIkhvbWUgUmFuZ2UgQW5hbHlzaXMgb2YgRWFzdGVybiBCb3ggVHVydGxlcyINCmF1dGhvcjogIjxjZW50ZXI+IE1vbGx5IFJpY2hhcmQgPC9jZW50ZXI+Ig0Kb3V0cHV0Og0KICBodG1sX2RvY3VtZW50Og0KICAgIGRmX3ByaW50OiBwYWdlZA0KICAgIHRvYzogeWVzDQogIGh0bWxfbm90ZWJvb2s6DQogICAgZGZfcHJpbnQ6IHBhZ2VkDQogICAgaGlnaGxpZ2h0OiBicmVlemVkYXJrDQogICAgbnVtYmVyX3NlY3Rpb25zOiBubw0KICAgIHJvd3MucHJpbnQ6IDEwDQogICAgdGhlbWU6IGNvc21vDQogICAgdG9jOiB5ZXMNCiAgICB0b2NfZmxvYXQ6DQogICAgICBjb2xsYXBzZWQ6IG5vDQogICAgICBzbW9vdGhfc2Nyb2xsOiB5ZXMNCiAgcGRmX2RvY3VtZW50OiBkZWZhdWx0DQplZGl0b3Jfb3B0aW9uczoNCiAgY2h1bmtfb3V0cHV0X3R5cGU6IGlubGluZQ0KDQotLS0NCg0KYGBge3Igc2V0dXAsIGluY2x1ZGU9RkFMU0V9DQprbml0cjo6b3B0c19jaHVuayRzZXQoZWNobyA9IFRSVUUpDQoNCnBhY2thZ2VzPC1jKCJhZGVoYWJpdGF0SFIiLCJkYXRhLnRhYmxlIiwiZ2dmb3J0aWZ5IiwiZ3JpZCIsIm1vdmUiLCJtb3ZlVmlzIiwiT3BlblN0cmVldE1hcCIsInBiYXBwbHkiLCJwbG90bHkiLCJyZ2RhbCIsInNwIiwidGlkeXZlcnNlIiwidmlyaWRpcyIpDQpzYXBwbHkocGFja2FnZXMsIHJlcXVpcmUsIGNoYXJhY3Rlci5vbmx5PVQpDQpgYGANCg0KDQpgYGB7ciBkYXRhLCBtZXNzYWdlPUZBTFNFLCB3YXJuaW5nPUZBTFNFfQ0KZGF0YSA8LSByZWFkX2NzdigiRWFzdGVybiBCb3ggVHVydGxlIFRyYWNraW5nLmNzdiIpDQoNCmRhdGEgPC0gZGF0YVshaXMubmEoZGF0YSQnbG9jYXRpb24tbG9uZycpLCBdDQpjb2xuYW1lcyhkYXRhKVtjb2xuYW1lcyhkYXRhKT09ImxvY2F0aW9uLWxvbmciXSA8LSAibG9uZyINCmNvbG5hbWVzKGRhdGEpW2NvbG5hbWVzKGRhdGEpPT0ibG9jYXRpb24tbGF0Il0gPC0gImxhdCINCmNvbG5hbWVzKGRhdGEpW2NvbG5hbWVzKGRhdGEpPT0iaW5kaXZpZHVhbC1sb2NhbC1pZGVudGlmaWVyIl0gPC0gImluZGl2aWR1YWwiDQpjb2xuYW1lcyhkYXRhKVtjb2xuYW1lcyhkYXRhKT09InV0bS1lYXN0aW5nIl0gPC0gImVhc3RpbmciDQpjb2xuYW1lcyhkYXRhKVtjb2xuYW1lcyhkYXRhKT09InV0bS1ub3J0aGluZyJdIDwtICJub3J0aGluZyINCmBgYA0KDQpFYXN0ZXJuIEJveCBUdXJ0bGUgTW92ZW1lbnQgaW4gTm9ydGggQ2Fyb2xpbmENCg0KRGF0YSBDb2xsZWN0ZWQgYnkgTi5DLiBXaWxkbGlmZSBDb21taXNzaW9uDQoNCkRhdGEgYWNjZXNzZWQgdGhyb3VnaCB3d3cubW92ZWJhbmsub3JnDQoNCg0KIyMjIFBsb3Qgb2YgUG9pbnRzDQpgYGB7ciBwbG90bHksIGVjaG89VFJVRSwgZmlnLmhlaWdodD02LCBmaWcud2lkdGg9OH0NCmdncGxvdChkYXRhPWRhdGEsIGFlcyh4ID0gbG9uZywgeT1sYXQsIGNvbG9yPWluZGl2aWR1YWwpKSsgDQogIGdlb21fcG9pbnQoKSsNCiAgICBsYWJzKHg9IkxvbmdpdHVkZSIsIHk9IkxhdGl0dWRlIikgKw0KICAgIGd1aWRlcyhjb2xvcj1ndWlkZV9sZWdlbmQoIklkZW50aWZpZXIiKSkNCmBgYA0KDQpgYGB7ciBsYXBwbHkgZnVuY3Rpb24sIG1lc3NhZ2U9RkFMU0UsIHdhcm5pbmc9RkFMU0UsIGVjaG89VFJVRSwgcmVzdWx0cz0naGlkZSd9DQpsYXBwbHkoc3BsaXQoZGF0YSwgZGF0YSRpbmRpdmlkdWFsKSwgDQogICAgICAgZnVuY3Rpb24oeCl3cml0ZS5jc3YoeCwgZmlsZSA9IHBhc3RlKHgkaW5kaXZpZHVhbFsxXSwiLmNzdiIpLCByb3cubmFtZXMgPSBGQUxTRSkpDQpgYGANCg0KYGBge3IgbGlzdCwgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRSwgZWNobz1UUlVFLCByZXN1bHRzPSdoaWRlJ30NCmZpbGVzIDwtIGMoIkJpbGx5IEJvYiAuY3N2IiwgIkJvbiBKb3ZpIChCSikgLmNzdiIsICJTQU1NWSAuY3N2IikNCmBgYA0KDQoNCmBgYHtyIGltYWdlcnksIG1lc3NhZ2U9RkFMU0UsIHdhcm5pbmc9RkFMU0UsIGVjaG89VFJVRSwgZmlnLmhlaWdodD02LCBmaWcud2lkdGg9OH0NCnV0bV9wb2ludHMgPC0gY2JpbmQoZGF0YSRlYXN0aW5nLCBkYXRhJG5vcnRoaW5nKQ0KdXRtX2xvY2F0aW9ucyA8LSBTcGF0aWFsUG9pbnRzKHV0bV9wb2ludHMsIA0KICAgICAgICAgICAgICAgICBwcm9qNHN0cmluZz1DUlMoIitwcm9qPXV0bSArem9uZT0xNyArZGF0dW09V0dTODQiKSkNCnByb2pfbGF0LmxvbiA8LSBhcy5kYXRhLmZyYW1lKHNwVHJhbnNmb3JtKA0KICAgICAgICAgICAgICAgIHV0bV9sb2NhdGlvbnMsIENSUygiK3Byb2o9bG9uZ2xhdCArZGF0dW09V0dTODQiKSkpDQpjb2xuYW1lcyhwcm9qX2xhdC5sb24pIDwtIGMoIngiLCJ5IikNCnJhc3RlciA8LSBvcGVubWFwKGMobWF4KHByb2pfbGF0LmxvbiR5KSswLjAwMSwgbWluKHByb2pfbGF0LmxvbiR4KS0wLjAwMSksIA0KICAgICAgICAgICAgICAgICAgYyhtaW4ocHJval9sYXQubG9uJHkpLTAuMDAxLCBtYXgocHJval9sYXQubG9uJHgpKzAuMDAxKSwgDQogICAgICAgICAgICAgICAgICB0eXBlID0gImJpbmciKQ0KDQpyYXN0ZXJfdXRtIDwtIG9wZW5wcm9qKHJhc3RlciwgDQogICAgICAgICAgICAgIHByb2plY3Rpb24gPSAiK3Byb2o9dXRtICt6b25lPTE3ICtlbGxwcz1XR1M4NCArdW5pdHM9bSArbm9fZGVmcyIpDQoNCmBgYA0KDQpgYGB7ciBpbWFnZXJ5IHBsb3QsIG1lc3NhZ2U9RkFMU0UsIHdhcm5pbmc9RkFMU0UsIGVjaG89VFJVRSwgZmlnLmhlaWdodD02LCBmaWcud2lkdGg9Nn0NCmF1dG9wbG90KHJhc3Rlcl91dG0sIGV4cGFuZCA9IFRSVUUpICsgdGhlbWVfYncoKSArDQogIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbj0iYm90dG9tIikgKw0KICB0aGVtZShwYW5lbC5ib3JkZXIgPSBlbGVtZW50X3JlY3QoY29sb3VyID0gImJsYWNrIiwgZmlsbD1OQSwgc2l6ZT0xKSkgKw0KICBnZW9tX3BvaW50KGRhdGE9ZGF0YSwgYWVzKGVhc3Rpbmcsbm9ydGhpbmcsDQogICAgICAgICAgICAgY29sb3I9aW5kaXZpZHVhbCksIHNpemUgPSAzLCBhbHBoYSA9IDAuOCkgKw0KICB0aGVtZShheGlzLnRpdGxlID0gZWxlbWVudF90ZXh0KGZhY2U9ImJvbGQiKSkgKyBsYWJzKHg9IkVhc3RpbmciLA0KICAgICAgICB5PSJOb3J0aGluZyIpICsgZ3VpZGVzKGNvbG9yPWd1aWRlX2xlZ2VuZCgiSWRlbnRpZmllciIpKQ0KYGBgDQojIyMgTWluaW11bSBDb252ZXggUG9seWdvbg0KYGBge3IgTUNQIHBsb3QsIG1lc3NhZ2U9RkFMU0UsIHdhcm5pbmc9RkFMU0UsIGVjaG89VFJVRSwgZmlnLmhlaWdodD02LCBmaWcud2lkdGg9Nn0NCm1jcF9yYXN0ZXIgPC0gZnVuY3Rpb24oZmlsZW5hbWUpew0KICBkYXRhIDwtIHJlYWQuY3N2KGZpbGUgPSBmaWxlbmFtZSkNCiAgeCA8LSBhcy5kYXRhLmZyYW1lKGRhdGEkZWFzdGluZykNCiAgeSA8LSBhcy5kYXRhLmZyYW1lKGRhdGEkbm9ydGhpbmcpDQogIHh5IDwtIGMoeCx5KQ0KICBkYXRhLnByb2ogPC0gU3BhdGlhbFBvaW50c0RhdGFGcmFtZSh4eSxkYXRhLCBwcm9qNHN0cmluZyA9IENSUygiK3Byb2o9dXRtICt6b25lPTE3ICtlbGxwcz1XR1M4NCArdW5pdHM9bSArbm9fZGVmcyIpKQ0KICB4eSA8LSBTcGF0aWFsUG9pbnRzKGRhdGEucHJvakBjb29yZHMpDQogIG1jcC5vdXQgPC0gbWNwKHh5LCBwZXJjZW50PTEwMCwgdW5vdXQ9ImhhIikNCiAgbWNwLnBvaW50cyA8LSBjYmluZCgoZGF0YS5mcmFtZSh4eSkpLGRhdGEkaW5kaXZpZHVhbCkNCiAgY29sbmFtZXMobWNwLnBvaW50cykgPC0gYygieCIsInkiLCAiaWRlbnRpZmllciIpDQogIG1jcC5wb2x5IDwtIGZvcnRpZnkobWNwLm91dCwgcmVnaW9uID0gImlkIikNCiAgdW5pdHMgPC0gZ3JpZC50ZXh0KHBhc3RlKHJvdW5kKG1jcC5vdXRAZGF0YSRhcmVhLDIpLCJoYSIpLCB4PTAuODUsICB5PTAuOTUsDQogICAgICAgICAgICAgICAgICAgICBncD1ncGFyKGZvbnRmYWNlPTQsIGNvbD0id2hpdGUiLCBjZXg9MC45KSwgZHJhdyA9IEZBTFNFKQ0KICBtY3AucGxvdCA8LSBhdXRvcGxvdChyYXN0ZXJfdXRtLCBleHBhbmQgPSBUUlVFKSArIHRoZW1lX2J3KCkgKyB0aGVtZShsZWdlbmQucG9zaXRpb249Im5vbmUiKSArDQogICAgdGhlbWUocGFuZWwuYm9yZGVyID0gZWxlbWVudF9yZWN0KGNvbG91ciA9ICJibGFjayIsIGZpbGw9TkEsIHNpemU9MSkpICsNCiAgICBnZW9tX3BvbHlnb24oZGF0YT1tY3Aub3V0LCBhZXMoeD1tY3AucG9seSRsb25nLCB5PW1jcC5wb2x5JGxhdCksIGFscGhhPTAuNSwgY29sb3IgPSAicmVkIiwgZmlsbCA9ICJyZWQiKSArDQogICAgZ2VvbV9wb2ludChkYXRhPW1jcC5wb2ludHMsIGFlcyh4PXgsIHk9eSkpICsgDQogICAgbGFicyh4PSJFYXN0aW5nIChtKSIsIHk9Ik5vcnRoaW5nIChtKSIsIHRpdGxlPW1jcC5wb2ludHMkaWRlbnRpZmllcikgKw0KICAgIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbj0ibm9uZSIsIHBsb3QudGl0bGUgPSBlbGVtZW50X3RleHQoZmFjZSA9ICJib2xkIiwgaGp1c3QgPSAwLjUpKSArIA0KICAgIGFubm90YXRpb25fY3VzdG9tKHVuaXRzKQ0KICBtY3AucGxvdH0NCg0KDQpwYmxhcHBseShmaWxlcywgbWNwX3Jhc3RlcikNCmBgYA0KDQoNCg0KIyMjIEtlcm5lbC1EZW5zaXR5IEVzdGltYXRpb24NCmBgYHtyIEtERSBwbG90LCBtZXNzYWdlPUZBTFNFLCB3YXJuaW5nPUZBTFNFLCBlY2hvPVRSVUUsIGZpZy5oZWlnaHQ9NiwgZmlnLndpZHRoPTZ9DQprZGVfcmFzdGVyIDwtIGZ1bmN0aW9uKGZpbGVuYW1lKXsNCiAgZGF0YSA8LSByZWFkLmNzdihmaWxlID0gZmlsZW5hbWUpDQogIHggPC0gYXMuZGF0YS5mcmFtZShkYXRhJGVhc3RpbmcpDQogIHkgPC0gYXMuZGF0YS5mcmFtZShkYXRhJG5vcnRoaW5nKQ0KICB4eSA8LSBjKHgseSkNCiAgZGF0YS5wcm9qIDwtIFNwYXRpYWxQb2ludHNEYXRhRnJhbWUoeHksZGF0YSwgcHJvajRzdHJpbmcgPSBDUlMoIitwcm9qPXV0bSArem9uZT0xNyArZWxscHM9V0dTODQgK3VuaXRzPW0gK25vX2RlZnMiKSkNCiAgeHkgPC0gU3BhdGlhbFBvaW50cyhkYXRhLnByb2pAY29vcmRzKQ0KICBrZGU8LWtlcm5lbFVEKHh5LCBoPSJocmVmIiwga2Vybj0iYml2bm9ybSIsIGdyaWQ9MTAwKQ0KICB2ZXIgPC0gZ2V0dmVydGljZXNocihrZGUsIDk1KQ0KICB2ZXIxIDwtIGdldHZlcnRpY2VzaHIoa2RlLCA3NSkNCiAgdmVyMiA8LSBnZXR2ZXJ0aWNlc2hyKGtkZSwgNTApDQogIGtkZS5wb2ludHMgPC0gY2JpbmQoKGRhdGEuZnJhbWUoZGF0YS5wcm9qQGNvb3JkcykpLGRhdGEkaW5kaXZpZHVhbCkNCiAgY29sbmFtZXMoa2RlLnBvaW50cykgPC0gYygieCIsInkiLCJpZGVudGlmaWVyIikNCiAga2RlLnBvbHkgPC0gZm9ydGlmeSh2ZXIsIHJlZ2lvbiA9ICJpZCIpDQogIGtkZS5wb2x5MTwtIGZvcnRpZnkodmVyMSwgcmVnaW9uID0gImlkIikNCiAga2RlLnBvbHkyPC0gZm9ydGlmeSh2ZXIyLCByZWdpb24gPSJpZCIpDQogIHVuaXRzIDwtIGdyaWQudGV4dChwYXN0ZShyb3VuZCh2ZXIkYXJlYSwyKSwiIGhhIiksIHg9MC44NSwgIHk9MC45NSwNCiAgICAgICAgICAgICAgICAgICAgIGdwPWdwYXIoZm9udGZhY2U9NCwgY29sPSJ3aGl0ZSIsIGNleD0wLjkpLCBkcmF3ID0gRkFMU0UpDQogIGtkZS5wbG90IDwtIGF1dG9wbG90KHJhc3Rlcl91dG0sIGV4cGFuZCA9IFRSVUUpICsgdGhlbWVfYncoKSArIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbj0ibm9uZSIpICsNCiAgICB0aGVtZShwYW5lbC5ib3JkZXIgPSBlbGVtZW50X3JlY3QoY29sb3VyID0gImJsYWNrIiwgZmlsbD1OQSwgc2l6ZT0xKSkgKw0KICAgIGdlb21fcG9seWdvbihkYXRhPWtkZS5wb2x5LCBhZXMoeD1rZGUucG9seSRsb25nLCB5PWtkZS5wb2x5JGxhdCksIGFscGhhPTAuNSwgY29sb3IgPSAicmVkIiwgZmlsbCA9ICJyZWQiKSArDQogICAgZ2VvbV9wb2x5Z29uKGRhdGE9a2RlLnBvbHkxLCBhZXMoeD1rZGUucG9seTEkbG9uZywgeT1rZGUucG9seTEkbGF0KSwgYWxwaGE9MC41LCBjb2xvciA9ICJncmVlbiIsIGZpbGwgPSAiZ3JlZW4iKSArDQogICAgZ2VvbV9wb2x5Z29uKGRhdGE9a2RlLnBvbHkyLCBhZXMoeD1rZGUucG9seTIkbG9uZywgeT1rZGUucG9seTIkbGF0KSwgYWxwaGE9MC41LCBjb2xvciA9ICJibHVlIiwgZmlsbCA9ICJibHVlIikgKw0KICAgIGdlb21fcG9pbnQoZGF0YT1rZGUucG9pbnRzLCBhZXMoeD14LCB5PXkpKSArDQogICAgbGFicyh4PSJFYXN0aW5nIChtKSIsIHk9Ik5vcnRoaW5nIChtKSIsIHRpdGxlPWtkZS5wb2ludHMkaWRlbnRpZmllcikgKw0KICAgIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbj0ibm9uZSIsIHBsb3QudGl0bGUgPSBlbGVtZW50X3RleHQoZmFjZSA9ICJib2xkIiwgaGp1c3QgPSAwLjUpKSArIA0KICAgIGFubm90YXRpb25fY3VzdG9tKHVuaXRzKQ0KICBrZGUucGxvdA0KfQ0KDQpwYmxhcHBseShmaWxlcywga2RlX3Jhc3RlcikNCmBgYA0KDQoNCg==